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ABSTRACT 

We have carried out a brown dwarf search over an area of 14' x 23' near the central por- 
tion of the Pleiades open cluster in five optical and near infrared bands {i' , Z, J, H, Ks) 
with 10 a detection limits of i' ^22.0, J ^20.0 and Ks ~18.5 mag. The surveyed area 
has large extinction in excess of Ay = 3 in the Pleiades region. We detected four 
new brown dwarf candidates from the colour-colour (J — K,i' — J) and the colour- 
magnitude (J, i' — K) diagrams. We estimated their masses as 0.046 Mq down to 
0.028 Mq. The least massive one is estimated to have a mass smaller than Roque 25 
ijMartin et al.lll998|) or int-pl-IZ-69 Ipobbie et al.ll2002b|) . and maybe the lowest mass 
object found so far in the Pleiades cluster. 
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1 INTRODUCTION 

Many brown dwarf (BD) searches have targeted young 
open clusters or star forming regions rather than general 
fields, because here BDs are still relatively bright and warm 
during their early phase. The Pleiades region is one of 
the best clusters for this purpose; it is a fairly young, 
nearby, rich and compact cluster. It ha s an age of 125 Myr 
iStauffer , Schultz, fe Kirkpatrick||l998l) at a distan ce of 130 
pc JCrawford fc Perrvl I197& IPinfield et alJl200oj) . Twelve 
hundred mem ber stars are located within 2.5 degre es of the 
cluster centre JPinfield. Jameson, fc Hodgkinlll99ai . 

A number of optical surveys searching for BDs in 
the Pleiades cluster have been done to date and iden- 
tified a numerous population of bona-fide BDs and BD 
candidates ; "CFHT-PL" objects found by an RI sur- 
vey dBouvier et al.l ll99Stk "NPL" objects by an RIJK 
survey (lFestidll99sf) ~ "MHObd" objects by a VI survey 
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"IPMBD" objects by an RI survey 
"Roqu e" objects by an IZ survey 
1999); "BPL" objects by an IZ 
, int-pl-IZ" objects by an IZ 
survey IIDobbi^et^l .I l2002bl) . About 20 of these objects are 
confirmed as Pleiades brown dwarfs on the basis of observed 
lithium abundance, radial velocity a n d prop er motion (e.g., 

119951 iRebolo et al 



iR^bclo^ajDatej^Osorip^&^Majrtinl 

19961: iBasri. Marcv. fc Grahamlll996l: I Zapatcro Osorio et al 



{)a 
~998t iMartm et al 

20001) . Ho wever, to dat e, onl; 
(Roque25; IMartm et al 



199 



* E-mail:chie@optik. mtk.nao.ac.jp; Present address: National 
Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, 
Japan 



one L-type BD candidate 
has been spectroscopically 
identified in the Pleiades cluster. Detection of L-type BDs 
is difficult at optical wavelengths because they should emit 
most of their bolometric flux in the infrared region. We 
have carried out a brown dwarf survey in both the optical 
and near infrared of a small region of the central region of 
the Pleiades cluster or southern portion of Merope. Most 
existing surveys have ignored this region because of its 
proximity to a small mo lecular cloud which results in large 
and variable extinction jBregerlll987l) . 
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Figure 1. The field locations on the sky. The shaded square 
shows the observed area and the open squares show the region 
covered only in optical observation. The overplotted dashed cir- 
cles are radii of 0.26° and 0.64° fro m the centre of the cluster, 
R.A.3 h 47 m ,Dec. +24°07' (J2000.0; iPinfield et ai1l200Cl) . 



2 OBSERVATIONS 

We surveyed an area of 14' x 23' centred at 
R.A.3 h 46 m 59 a .6,Dec. + 23°44'59"(J2000.0) both in the 
optical and near infrared bands. The observed area is shown 
as the shaded rectangle in Fig. The four open rectangles 
show the area covered at optical wavelengths. 



2.1 Near Infrared Survey 

Our near infrared survey was made on 2000 October 12, 17 
and 18 with SIRIUS (Simultaneous InfraRed Imager for Un- 
biased Survey) on the University of Hawaii 2.2 m telescope, 
atop Mauna Kea. SIRIUS employs three 1024 x 1024 pixel 
HgCdTe arrays which provides J (1.25^m), H (1.65/xm), and 
Ks ( 2 .15Atm) band images sim ultaneously jNagashima et all 
119991: iNaeavama et alJ l2002). The field of view in the each 
band is 4'9 x 4!9 with a pixel scale of C// 28 at the Cassegrain 
focus of f/10. We mapped an area of 14' x 23' with 3x5 
tiles of the field of view. We obtained each of the tiles by 
dithering 18 frames. The exposure time was 60 sec for each 
frame. The limiting magnitudes at 10 a are 20.0, 19.2, and 
18.5 at J, H, and Ks, respectively. Typical seeing in the Ks 
band was l'.'O. 

Dark and dome-flat frames were taken at the beginning 
or the end of each night. We observed standard stars, Nos. 
9105, 9107, 91 08, 9116 and 9 1 88, in the faint NIR standard 
star catalog of IPersson et alJ il998f) on the same nights for 
photometric calibration, and the red standard stars, LDN 
547, BRI B0021-02 14 and BRI B2202-1119, in Table 3 of 
IPersson et alJ {l998) on 2001 August 31 for determination of 
the colour conversion formula from ours to the CIT system. 
We obtained the following equations as colour-conversion 
from ours to the CIT system. 

Jcit = Jobs-(0.013±0.017)x(J-A"s) obs + (0.012±0.030) 
#ciT = #ob s -(0.029±0.029)x(J-#) obs + (0.007±0.033) 
7vciT = A"s obs -(0.024±0.020)x(J-A"s) obs + (0.018±0.035) 

We used pipeline software based on NO AO's IRAF 



(Imaging Reduction & Analysis Facility) 1 package to re- 
duce the data. We have applied standard NIR image re- 
duction procedures, including dark current subtraction, sky 
subtraction and flat fielding. Source detection, morphologi- 
cal classification and photometry were performed with the 
APPHOT package in IRAF. 

2.2 Optical Survey 

The i' (0.77 fim, SDSS system) and Z (0.9 fim, RGO) 
band image was taken on 2000 October 20 with the Wide 
Field Camera (WFC) on the Isaac Newton Telescope (INT) 
at La Palma, as par t of the INT Wide Angle Survey 
jMcMahon et al.ll200ll) . The WFC accommodates four 2048 
x 4196 pixel EEV CCDs. The field of view for each CCD is 
ll!4 x 22'8 with a pixel scale of 0'.'33 at the prime focus of 
the INT and the mosaic layout of the four CCDs enables us 
to encompass the whole field of 14' x 23' taken in the near 
infrared bands (Fig. . We obtained one frame for each fil- 
ter; the exposure time of each frame was 600 s. The 10 a 
detection limits are 22.0 and 20.5 at i 1 and Z, respectively, 
where i! and Z are on natural system of WFC. 

The observed data were reduced at the Institute of 
Astronomy, Cam bridge, using the W FC data reduction 
pipeline software jlrwin fc LewislfeoOlf) . All the frames were 
bias-subtracted and then corrected for linearity before the 
science fields were flatfielded and defringed. The pipeline 
was also used to generate photometrically and astrometri- 
cally calibrated object catalogues. On the night in question, 
the %'- and Z-band photometric zeropoints were measured 
from three s tandard s tar fi elds (observed once each) in the 
catalogue of lLandoit] J1992T) and were found to be stable to 
about 1 %. The astrometric solution takes account of the 
radial distortion across the WFC field of view, and resulted 
in residuals of 0'.'22. 



3 SELECTION OF PLEIADES BD 
CANDIDATES 

3.1 The colour-colour diagram 

We constructed the colour-colour (J — K, i' — J) diagram for 
all the objects, 846 sources in total, detected with S/N ratios 
higher than 10 in the near infrared and optical bands, dif- 
ferentiating them as extended (149; open circles) and point 
(697; filled circles) sources in Fig. |U Seven sources detected 
with S/N ratios higher than 10 in the near infrared but 
lower than 5 in the optical are also plotted with large error 
bars. "Extended objects" were sorted with a criterion of the 
FWHMs larger by at least 20 % than the seeing size of the 
night in the Ks images. 

The thick solid and dashed curves from the bottom-left 
corne r represents the loc i of main sequence and red giant 
stars tBe ssellll99ll . ll979]) : spectral types of K7 dwarf at the 
cusp at around (1,0.8) and of M7 at the tip (3.5,1.2) for (if — 
J, J — K) , while the arrow at the upper-right corner shows 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 2. The J — Ks, i' — J colour-colour diagram for all objects detected in all the 5 bands. Small dots are point sources and open 
circles are extended sources. Asterisks are the known late-type field dwarfs. The thick solid and dashed curves from the bottom left 
corner represent the loci of main sequence and red giant stars. The arrow at the upper right corner shows the reddening vector of Ay = 
3 mag. The thin dashed line is the ~125 M yr isochrone of the N extGen model fearaffe et aljll99St) . and the thin dot-dashed line is the 
~120 Myr isochrone of the Dusty model iChabrier et alj|2000l) . Numbers with tick marks denote masses in solar mass units for each 
model. 



the reddening vector of Ay = 3 mag ijSchultz fc Wiemerl 
Il975l: iRieke fc Lebofskvlll985l) . Most of the point sources 
(filled circles) are concentrated above the locus from K7 to 
M7, along the extinction vector up to 3 mag from the main 
sequence (reddening band) , indicating they are K- or M-type 
dwarfs lying intrinsically on the main sequence. On the other 
hand, most of the extended sources and some point sources 
spread upward of the reddening band. About 20 of the ex- 
tended sources have galaxy-like shapes. The K-correction 
curves from z = - 0.1 of elliptical and spiral galaxies (Fu- 
rusawa, private communication) are plotted as dotted lines. 
Thus, objects spread upward of the reddening band includ- 



ing the point sources (182 objects in total) are likely to be 
elliptical or spiral type galaxies reddened with K-corrections 
of z ~ to 0.1 and Ay of to several magnitudes. A major 
problem in this kind of survey is to dist inguish stars from 
distant galaxies at faint magnitudes (e.g. lFestinlll997r) . The 
optical and near-infrared colour-colour diagram allows us to 
easily distinguish them. Note that the near infrared colour- 
colour diagram (J — H,H — K) alone cannot distinguish 
low-mass stars from t hese galaxies, due to merging of them 
in th e diagram (e.g. iKirkpatrick et al"lll999l : IJarrett et alJ 
l2000h . 

To highlight the expected location of Pleiades BDs in 
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the colour-colour diagram, two theoretical isochrones for 
the Pleiades cluster are also plotted in Fig. [5] the dashed 
line is the 125 My r isochrone of the NextGen model 
iBaraffe et alj fl998"l. while the dot-d ashed line is the 120 
Myr isochrone of the Dusty model JChabrier et al.ll200l 
; Baraffe, private communication for 120 Myr model). The 
7c magnitudes of these models have been transformed on to 
our system using the following relation, i' = 7c + 0.211(7?c — 
7 C ) + 0.011(V - 7 C ) 2 - 0.003(1/ - 7 C ). This transforma- 
tion is derived via Landolt system (7L an ); the transfor- 
mation between i' and Ihun is from a n INT web p age 2 , 
and that between 7L an and 7c is from Bessell (1990). The 
NextG en model is based on the nongray dust-free atmo- 
sphere llHauschildt. Allard. fc Baronl I1999T) , and describes 
successfully various observed properties of early-mid M 
dwarfs, but predicts near infrared colours too blue compared 
to the observed ones for cooler objects. The Dusty model 
includes the effect of grain formation, that is, 1) the pho- 
tospheric depletion of dust-forming elements, and 2) scat- 
tering and absorption by dust, explaining successfully the 
near infrared colours of the late M and L t ypes with 1,800 
< T off < 2, 200 K fe.g. lJameson et al.lEoO^t . Because of the 
inadequacy of the NextGen model in the cooler BD region, 
we have adopted the Dusty model in the following discus- 
sion. Using near infrared p hotometry and lithium tests o f 
low mass cluster members, IZapatero Osorio et al] Jl997bT) 
estimated the stellar / substellar boundary to be located at 
an absolute magnitude of (~ 0.075 Mq) for Pleiades (~120 
Myr) to be at Mi = 12.4, Mj = 10.0 and M K = 9.0 mag. 
This corresponds to the point (i — J, J — K) = (2.9, 1.0) 
in our colour-colour diagram, highlighted by a cross, at ap- 
prox imately spectral type M6. Twelve late-type field dwarfs 
from lDobbie et all (|2002ar ) were also plotted as asterisks. We 
adopt the redder side of the line running from the cross in 
parallel with the extinction vector as the location of Pleiades 
BDs in Fig. |21 Note that this location coincides with that of 
late-M and L-type field dwarfs. There remain 12 sources on 
the redder side of the line. 



3.2 The colour-magnitude diagram 

To separate Pleiades BDs from field dwarfs, we construct 
the colour-magnitude diagram (J, %' — TsT). The 12 red ob- 
jects selected by Fig.|5|are plotted in Fig. [31 The theoretical 
tracks representative of the low mass stellar and substellar 
members of the Plei ades are also plotted; N extGen model 
fearaffe et al] 1 19981) and Dusty model iChabrier et al] 
2000). We choose the J magnitude as the ordinate of the di- 
agram because i' magnitudes of the red objects have larger 
errors and K magnitudes are less sensitive to the mass differ- 
ence on the theoretical tracks. The arrow at the upper-right 
corner indicates the reddening vector of Av = 3 mag, which 
runs nearly in parallel with the tracks. 

We adopt the distance modulus used bv lDobbie et al] 
(2002b), (m— A7)o = 5.53, to estimate the ordinate displace- 
ment (J magnitude) of the isochrone. Different age estimates 
(70 - 150 Myr) result in a displacement of ~ -0.3 and ~ +0.1 
mag. Different distance estimates result in a displacement of 
~ -0.2 mag (Hipparcos) and ~ +0.1 mag (photometry). The 

2 http:/ /www. ast.cam.ac.uk/~wfcsur/photom.php 



effect of the cluster depth results in a displacement of ~ ± 
0.2 mag. Additionally, location of unresolved binaries is by 
0.75 mag above the single-star sequence (e.g. iPinfield et al] 
2000). Accounting for all of these uncertainties and allowing 
for a small degree of error in both of the theoretical mod- 
els and photometry, BDs belonging to the Pleiades cluster 
should be located between 0.3 mag below and 1.0 mag above 
the relevant single-star isochrone. 

With these criteria, we identified four Pleiades BD can- 
didates from 12 red objects. The remaining eight objects 
lie more than 1 mag below the isochrone in the colour- 
magnitude diagram, although they fall in the Pleiades BD 
domain in the colour-colour diagram. They could be late-M 
and L-type dwarfs lying background of Pleiades. Finally, we 
list the four sources (PL- # of Fig.[3J as BD candidates from 
our five-colour survey. All of them are previously unpub- 
lished. The positions and magnitudes of four BD candidates 
are provided in Table and that of the remaining eight ob- 
jects are in Table |5] Finding charts of four BD candidates 
are presented in Fig. 2] 



4 DISCUSSION 

4.1 Estimation of the BD candidate masses 

To estimate masses of BD candidates from their observed 
magnitudes, we should determine the level of interstellar ex- 
tinction towards them. The distribution of the background 
field stars in Fig. [5] indicates that the interstellar extinc- 
tion in the surveyed portion of Pleiades ranges from to 
3 mag in Av- The foreground extinction toward Pleiades is 
relatively uniform with the value of Ay ~ 0.12 mag, except 
for a region south of Merope, which contains a small molec- 
ular cloud llCrawford fc Perrvlll976l) . Because the cloud is 
assigned to be inside the cluster from polari metry and red - 
dening of the member stars and the field stars ]Bregerll987fl . 
the extinctions estimated from background stars should be 
regarded as upper limit for the BD candidates. 

To estimate the maximum interstellar extinction for the 
individual BD candidates, we constructed an extinction map 
in the surveyed area from all the field stars (small dots be- 
tween the two straight lines in Fig. |5J , assuming them to be 
K- or M-type main sequence stars (the thick line in Fig.|2J. 
Visual extinction, Av, toward each star was calculated from 
its i' — J and J — K c olour s with the reddening la w of 
ISchultz fc Wiemerl <1975h and lRieke fc Lebofskvl (1198511 . i.e. 
E(J - K)/E(i' — J) = 0.56, A v = 6.10 x [(J - K) obs — (J — 
if) intrinsic]. We meshed the area with 17 x 10 grids of 1' 
square and adopted the maximum Av within the grid, and 
averaged them over each mesh with Gaussian weights. The 
final Ay map is shown in Fig. |H| and the values toward the 
BD candidates are listed in Table |3] 

The estimated maximum extinctions for individual can- 
didates are shown in Fig.[3]with dotted arrows. If extinction 
is neglected, PL-2 and -3 a re estimated to hav e M ~ 0.035 
Mq, similar to Rogue 25 ([M jxtrn et alj |l998T) and int-pl- 
IZ-69, -81, -84 jDobbie et al.ll2002bl) . the coolest candidate 
members so far detected. On the other hand, assuming the 
maximum likely level of extinction, their masses are esti- 
mated to be M ~ 0.04 Mq. Similarly we estimate the mass 
of the faintest candidate to lie in the range of 0.033 - 0.028 
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Figure 3. The J,i' — K colour-magnitude diagram. The BD candidates are shown with filled circles. Dotted arrow from each candidates 
shows the estimate d maximum reddening vector (see Section 14.11 . Other symbols are as in Fig. El The i' - K colour of Roque 25 
jMartin et a l. 1998) is also shown with a perpendicular line due to lack of J magnitude. 



Table 1. The positions and magnitudes of new brown dwarf candidates. 



ID 


R.A. 


Dec. 


i' 


Z 


J 


H 


K 


name 


c06 
c08 
clO 
cl2 


03:45:50.6 
03:46:34.3 
03:45:11.7 
03:45:58.5 


+23:44:37 
+23:50:04 
+23:41:44 
+23:41:54 


20.02 ± 0.02 
21.36 ± 0.06 
21.98 ± 0.09 


18.19 ± 0.01 
19.52 ± 0.03 
19.88 ± 0.05 
21 69 + 029 


16.34 ± 0.06 

17.35 ± 0.06 
17.59 ± 0.06 
18.54 ± 0.08 


15.68 ± 0.07 
16.63 ± 0.07 
16.79 ± 0.07 
17.42 ± 0.09 


15.19 ± 0.07 
16.01 ± 0.08 
16.13 ± 0.08 
16.53 ± 0.09 


PL-1 
PL-2 
PL-3 
PL-4 



Table 2. The positions and magnitudes of other red objects. 



ID 


R.A. 


Dec. 


i' 


Z 


J 


H 


K 


cOl 


03:46:24.9 


+23:46:10 


22.28 t o.ll 


20.93 ± °il 


19.30 ± 0.08 


18.90 ± 0.12 


18.28 ± 0.13 


c()2 


03:45:30.9 


+23:49:51 


22.03 ± 0.1 


20.45 ± 0.08 


18.97 ± 0.07 


18.32 ± 0.08 


17.91 ± 0.09 


c03 


03:45:59.7 


+23:40:53 


22.33 ± 0.14 


20.82 ± 0.12 


18.93 ± 0.07 


18.35 ± 0.08 


17.89 ± 0.10 


c04 


03:46:26.6 


+23:39:35 


21.57 ± 0.07 


20.14 ± 0.06 


18.01 ± 0.06 


17.22 ± 0.07 


16.72 ± 0.08 


c05 


03:46:28.0 


+23:42:39 


22.19 ± 0.12 


21.11 ± 0.15 


18.60 ± 0.07 


17.74 ± 0.08 


17.29 ± 0.09 


c07 
c09 
ell 


03:46:20.7 
03:45:50.9 
03:45:49.2 


+23:41:56 
+23:44:57 
+23:41:26 


23.18 + lf 7 
22.65 t°o.\l 
24.34 tli 


22.18 + °;<g 
20.84 + °™ 
21 Q3 + 026 


19.14 ± 0.08 
18.63 ± 0.07 
19.88 ± 0.10 


18.75 ± 0.10 
17.64 ± 0.08 
18.78 ± 0.13 


18.03 ± 0.11 
17.08 ± 0.09 
18.33 ± 0.14 
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Figure 4. Finding charts for new brown dwarf candidates (Ks 
band). Each panel is 2' X 2', North is up, East is left. 




An [magi 

Figure 5. The extinction map constructed from the J — K and 
i' — J colours of the point sources in the reddening band. Small 
dots are field stars and diamonds are BD candidates. Estimation 
of the extinction is described in Section 14.11 The correspondence 
between grey scale and Ay is indicated in the scale bar. 



Mq. Therefore, PL-4 maybe the lowest mass member of the 
Pleiades identified to date. The estimated masses from J 
magnitudes for the case of no extinction and that of max- 
imum possible extinction toward individual candidates are 
shown in Table [3] Note that uncertainties in the age and the 
distance of the cluster and the effect of the cluster depth 
contribute an additional error of 10 - 20 % in the mass es- 
timation; e.g. the mass of PL-4 for the case of no extinction 
is 0.0281 'qqj Mq. We note that our mass estimates would 
increase slightl y if we were to take in to account the "missing 
M dwarf gap" foobbie et al.ll2002rj) . 



Table 3. Mass estimates for brown dwarf candidates based on 
120 Myr DUSTY model. 







Mass(MQ) 




Name 


MAX A v 


A v =0 - MAX 


SpT 


PL-1 


2.5 


0.046 - 0.063 


M6 - M8 


PL-2 


2.0 


0.035 - 0.040 


M8 - LI 


PL-3 


2.5 


0.033 - 0.039 


M8 - LI 


PL-4 


3.5 


0.028 - 0.033 


LI - L5 
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Figure 6. The broad-band energy distributions for four brown 
dwarf candidates. Solid and dashed lines represent brown dwarf 
candidates with no extinction and with possible maximum ex- 
tinction, respectively. M6 - M8 and L0 - L5 field dwarfs are over- 
plotted with dotted lines. 



4.2 Broad-band energy distributions 

In Fig. |S| we compare the four SEDs with those of cool 
field dwarfs, determined in previous surveys (M6 - M8; 
Kirkpatrick & McCarthy 1994, L0 - L5; Kirkpatrick et al. 
1999,2000). Solid and dashed lines represent SEDs of the 
four BD candidates for the cases of no extinction and max- 
imum extinction discussed in § 14.11 respectively. We have 
assu med the distance modulu s for Pleiades of (m — M)o = 
5.53 (ICrawford k, Perrvlll97ri iPinfield et al"ll2000D to con- 
vert the observed magnitudes to absolute magnitudes for 
the cluster member. The zero-mag nitude fluxes are taken 
from lBessell. Castelli. fc Pled Jl998i) . 

All the candidates show a steep decline towards the 
shorter wavelengths and a broad hump between 1 and 2 
/jm. Within the uncertainties arising from the correction for 
interstellar extinction, the four SEDs of the candidate BDs 
are well matched by the energy distributions of the late-M 
and L-type field dwarfs: PL-1 is consistent with a mid to 
late-M spectral type, PL-2 and PL-3 with late-M to early L 
spectral type and PL-4, an early-mid L spectral type. We 
summarize these results in Table |3" 
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4.3 Contamination of field stars 

As discussed in section 13.11 the optical and near-infrared 
colour-colour diagram allows us to separate distant galaxies 
and giants from BD candidates. However, field dwarfs are 
still possible sources of contamination. 

The Dusty model indicates that Pleiades BDs are over- 
luminous with respect to field stars of similar spectral type 
by ~ 1 mag. The criterion of selecting objects between 0.3 
mag below and 1.0 mag above the isochrone means that we 
are sensitive to single and binary dwarfs at distance of 71 
- 130 and 96 - 184 pc, respectively. This corresponds to 
space volumes of 15.5 and 45.1 pc 3 for single and binary 
field stars, respectively. Ass uming a space de nsity of M8 - 
L4.5 dwarfs of 0.0066 pc~ 3 JGizis et alJ|2000h and a binary 
fraction of 50 % (e.g. ISteele fc Jamesonlll995l) . we estimate 
the contamination of ~ 0.2 field star in our BD candidates. 
Thus, the level of the contamination of the BD candidates 
by non-members is negligibly low. 



n{n) 



= 12.78. 



(2) 



iVtotal 

iVannul [n(r ) - n(r z )] 

Finally, we calculate the values of substellar mass func- 
tion of the cluster, iVtotal per unit mass. The values are 
given in Table HI a nd are plotted in F i g. |7| with several pre- 
vious estimates dDobbie et al J l2002bt iHodekin fc Jameson! 
l2000t lHamblv et al.lll999ft . The shape of mass function of 
the Pleiades across and below the stellar / substellar bound- 
ary has been reported to be well represented by a power- 



investigator (e.g. Jameson ct al. 2002; 


Bouvier et alJll998 


IHodekin & Jameson! 


2000 




ZaDatero Osorio et alj 


1997a) 



Our new points seem to support the larger values of a. How- 
ever, as the present survey covers only a small fraction of 
the cluster (i.e., < 0.5 %), this result must be tr eated with a 
measu re of caution. Nevertheless, we note that iBeiar et alJ 
d200ll) have recently determined the substellar mass function 
of the much younger a Orionis cluster to be well matched 
by a power-law with index a = 0.8. 



4.4 The cluster luminosity function and mass 
function 

We have derived the cluster luminosity function for J magni- 
t udes between 15.5 to 20.0 which correspond to the BD limit 
o flZapatero Osorio et al.l dl997bT) and the 90 % completeness 
limit of our survey (almost the same as 10 a detection limit), 
respectively. In our calculation we assume that the interstel- 
lar extinction toward Pleiades in our survey area is spatially 
uniform with the value of 1.5 mag at V or 0.4 mag at J, we 
subdivided this magnitude range into 2-mag bins, given in 
Table H 

Using the Dusty model, we have calculated the masses 
corresponding to each of our luminosity bins. We estimated 
the total number, iV annu i, of member of the whole clus- 
ter from the number of our sample. First we calculate the 
number, iVannui, of members inside an annulus. Our survey 
area covers 7.7 % of the area of the annulus of an inner 
and an outer radius of r, = 0.26° and r — 0.64°, respec- 
tively (Fig0. 7V an nui is calc ulated simply from the area ra- 
tio. N ext we calculate iVto tai ■ |PmfiekLJa meson. fc Hodekinl 
(1998) and lRaboud fc Mermilliodl (Il998l) show that the sur- 
face density distribution of the Pl eiades mem bers can be 
well fitted by a King distribution (lKing|ll"962l) whose core 
radius increases as the stellar mass decreases. The following 
equation, which is obtained by integrating the function of 
King distribution, provides the number of cluster members 
inside the circle of a radius r, 



n(r) 



ln(l + a;)-4 



vO+jb- 1 



+ 



l + x t 



(1) 



where x = (r/r c ) 2 , Xt = (rt/r c ) 2 , k is the surface density of 
the cluster, r the radius from the cluster centre, r c the core 
radius at which the surface density of members falls to half 
its central value, and rt the tidal radius at which the exter- 
nal gravitational potential equals the cluster potential. The 
total number of cluster members, iVtotal > can be determined 
independently of k by taking the ratio n(rt)/[n(r ) — n(ri)], 
where n(r ) — n(r;) is equal to iVannui and n (r t ) is equal to 
iVtotal - We adopt a value of r c — 2.2° from Ijameson et alJ 
(2002) for low- mass members of Pleiades, and a v alue of 
r t = 5.54° from lPinfield. Jameson, fc Hodekinl Em 



5 CONCLUSIONS 

We have undertaken an optical (i',Z) and near infrared 
(J, H, Ks) survey of a 14' x 23' area near the centre of the 
Pleiades cluster to search for brown dwarf members. We 
have unearthed 4 BD candidates (PL-1, -2, -3, and -4), out 
of 853 sources. PL-1 is likely an M-type BD, PL-2 and -3 are 
likely late-M or early-L type BDs, while PL-4 is probably an 
L-type BD. The estimated mass of our faintest candidate, 
PL-4, is in the range of 0.028 - 0.033 M , lower than that of 
Roque 25 dMartin et aljfl 998) and can be lower than that of 
int-pl-IZ-69 dDobbie et al.l2002bf) . possibly making it lowest 
mass Pleiad identified to date. 
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